Backreaction of accreting matter onto a black hole in the 
Eddington-Finkelstein coordinates 

E. Babichev a,fe , V. Dokuchaev 6 , Yu. Eroshenko b 

a Laboratoire de Physique Theorique d'Orsay CNRS UMR 8627, 
■ Universite Paris-Sud 11 91405 Orsay Cedex, France, 

b Institute for Nuclear Research of the Russian Academy of Sciences, 
60th October Anniversary Prospect, la, 117312 Moscow, Russia 



Abstract 



o 

<N 
X> 

• We study backreaction of accreting matter onto a spherically symmetric black hole 

in a perturbative way, when accretion is quasi steady-state. General expressions for 
corrections to the metric coefficients are found in the Eddington-Finkelstein coordi- 
| nates. It is shown that near the horizon of a black hole, independently of the form of 

the energy-momentum tensor, the leading corrections to the metric are of the Vaidya 
| ' form. Relation to other solutions is discussed and particular examples are presented. 

Mi 

1 Introduction 

>. 

. Study of spherically-symmetric black hole accretion has a long story. In Newtonian gravity, 

first the problem of matter accretion onto compact objects has been formulated in a self- 
■ similar manner by Bondi [lj. Later, Michel investigated steady-state spherically symmetric 

. flow of a test gas onto a Schwarzschild black hole in the framework of General Relativity 

[2]. Accretion of a perfect fluid with a general equation of state onto a Schwarzschild black 
hole has been investigated in [31 SJ [5], and a similar analysis for a charged black hole has 
been done in [6] . In [7] an exact analytic solution of non-spherical accretion of a stiff fluid 
onto a rotating black hole was derived, and a generalization of this analysis to the case 
^ " of a rotating charged black hole was done in [8j. Other studies of spherically symmetric 

accretion of different types of fluids onto black holes and wormholes were undertaken in 
a number of works \TU[ \TT\ I12j . See also an excellent review on black holes in FRW 
universe [13]. 

Usually in analytic study of the problem of accretion an important simplification is 
made — backreaction of an accreting matter is neglected. In practice this assumption 
means that the accreting matter is very light. It is, however, of a great interest to inves- 
tigate backreaction effects on the metric of a black hole. To solve this problem in a full 
way is notoriously difficult, only a few analytic solutions are known, where backreaction 
is taken into account. In particular, the celebrated Lemaitre-Tolman solution describes a 
spherically symmetric collapsing dust [141 [15] . Another well-known solution is the Vaidya 
solution [16\ [T71 [TBI [19] , which can be applied for description of a black hole surrounded 
by radially moving infalling photons. Generalizations of the Vaidya solution to include 
more general energy-momentum tensor are possible, see [20J. 

In this paper we take a different approach to attack the problem of backreaction of 
accreting matter. Instead of looking for particular exact solutions, we use a perturbative 
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scheme: for an accreting matter with general energy-momentum tensor we find corrections 
to the metric due to backreaction. The scheme goes as follows. Zero approximation is a 
solution for an accreting matter in the test fluid approximation, i.e. when backreaction 
of accreting matter is neglected. The first correction to the metric are found from the 
Einstein equations, where the source term is computed using zeroth-order solution for 
matter. Although the perturbative approach is not designed to find exact solutions, neither 
to take into account strong backreaction effects, it allows to study the problem in general 
manner, without relying on a specific choice of the energy-momentum tensor. We also 
note, that a similar approach was implemented in [21] for calculating backreaction effects 
using the Schwarzschild (static) coordinates. Here, instead, we use ingoing Eddington- 
Finkelstein coordinates, in which case we avoid problems related to the singularity of the 
coordinate system on the horizon of a black hole. 

The paper is organized as follows. In Section [2] the general perturbation scheme is 
described, Sec. 12. 1| then, according to this scheme, corrections to the metric are found 
for the case of accretion onto a Schwarzschild black hole, Sec. 12. 2\ and onto a Reissner- 
Nordstrom black hole, Sec. 12.31 The results are applied to accretion of a perfect fluid with 
an arbitrary equation of state, Sec. 12.41 The shift of the apparent horizon is discussed 
in Sec. 12.51 In Section [3] we compare corrections found with the help of the perturbative 
approach to known exact solutions. It is emphasized that in the near horizon limit, 
corrections to the metric are of the Vaidya form, independently of the energy-momentum 
tensor of an accreting matter. Section H] contains particular examples: corrections to the 
metric are computed for accretion of a stiff fluid, Sec. I4.lt and for accretion of a canonical 
scalar field and a galileon field, Sec. 14.21 We conclude in Sec. [5j 

2 Corrections to the metric 
2.1 Perturbation scheme 

In general case of a black hole accretion the full solution to the Einstein equations are not 
known. Indeed, there are only a few special cases of the energy-momentum tensor, for 
which the exact solutions were found, e. g., Tolman [15] and Vaidya |16| solutions. Here 
we proceed with the perturbative approach to the problem. We find corrections to the 
metric due to accretion of matter with a general energy-momentum tensor, Tau, in the 
(almost) steady-state regime of accretion. In this regime the energy-momentum tensor is 
a function of the radial coordinate r only. The solution for the energy-momentum tensor 
Tuu(f) is assumed to be found in the test fluid approximation, i.e. without taking into 
account backreaction. Then, using T„ v {r) as the zeroth-order approximation, we can find 
the first order correction to the metric due to the back reaction of the energy-momentum 
tensor. 

More precisely, our approach is as follows. The full set of equations consist of the 
Einstein equations, 



where by <p we denote collectively the degrees of freedom for the accreting matter. For 
example, in the case of accretion of a single scalar field, <j) is a scalar field itself. In fact, the 




(1) 



and the equations of motion for matter 
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equations of motion for matter ([2]) can be obtained from the Bianchi identities, therefore 
they are not independent, but in practice it is convenient to keep them. 

A zeroth-order approximation is found by neglecting backreaction of the matter in 
Eq. ([1]). Put differently, the solution in this approximation is the vacuum solution for the 
g^f, so that 



metric, g$ 



G 



0. 



(3) 



In the same approximation, the solution for the matter field(s) 



is found from, 



£ 



„(o) 



0. 



(4) 



Solutions in this zeroth-order approximation has been used in literature extensively. 

A step beyond the test-fluid approximation is to find corrections to the metric due to 
backreaction of matter. In a first-order approximation we can substitute <p^ and gff} in 
the r.h.s. of (pQ), so that the Einstein equations are of the form, 



G 



fJ,U 



JO) + (1) 



8itT, 



fJIV 



(5) 



Assuming that gfy is small with respect to g^J one can linearize the above equation to 
obtain a system of linear partial differential equations. We will implement this approach 
below. 



jo) 



2.2 Schwarzschild black hole 

The general spherically symmetric line element can be written in the following form, 
containing two arbitrary functions v(V,r) and X(V,r) (see, e.g., [22l [23] ) : 

ds 2 = e »(V,r) + 2X(V,r) dv 2 _ 2e KV,r) d y dr _ ^ (g) 

This is a coordinate frame of the Eddington-Finkelstein (EF) type, related to the infalling 
null (photon) geodesies. Note that the metric Q is similar to the one introduced by H. 
Bondi et. al. [23] using "radiation coordinates". The vacuum Schwarzschild solution is 
recovered by setting A = and e u ^ V ' r ^ = 1 — 2Mq/t, where Mq is the mass of the black 
hole, so that 

dg2 = / _ 2A^ dy2 _ _ (7) 

Note that in our conventions Eq.([7]) is the solution for the metric in the zeroth-order 
approximation. 

In analogy to the line element (J7]), instead of the metric coefficient i/(V,r) we will use 
for convenience a function M(V,r), defined as follows, 

e*™ = 1 - (8) 
r 

so that in the zeroth-order approximation, M(V, r) = Mq = const. 
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Plugging ([8]) and (|6|) into the perturbed Einstein equations © we obtain the following 
system of equations (other components of © are trivial), 

8^0° = -^ + 7)+^ 0) 

SttTq 1 = yZ>, (10) 
r 

8^ = -e^^ + -V^--^, (12) 



SkT 2 2 = 8vrr 3 3 = -e v ( A" + — ) - e~ x X' 

-eMA« + £ + ^ + W|. (13) 



where dots stand for d/dV, and primes denote d/dr. Note that according to the general 
scheme, see Eq. ([5]), the l.h.s. of Eqs. ([9lll3p contain components of the energy- momentum 
tensor, taken in the zeroth-order approximation, i.e. on solutions of the matter field(s) 
when backreaction is neglected. Not all components of the above system of equation are 
independent — using the Bianchi identities it can be shown that Eq. (|13p is a combination 
of Eqs. (f9j)- (fT2jh Thus we will not take Eq. (fl~3j) into account below. 
Substituting © into (J9j) and jTO]) we find, 



M / = 47rT V (14) 

and 

M = A, (15) 
correspondingly, where in the last equation we introduced a notation 

A=-A^T l r 2 , (16) 

for the total flux of energy crossing some radius r. The r.h.s. of Eqs. (|14p and (|15|) 
are taken in the zeroth-order approximation. In this approximation, in particular, T Q ° ~ 
Tq\ and Tq 1 ~ Too- Also, the components of the energy-momentum tensor are time- 
independent in the given approximation. The flux A does not depend on r as well, i.e. 
A =const. For each particular case of black hole accretion the value of A is calculated 
in the zeroth-order approximation by solving the equations of motion in a fixed metric 
(normally by finding integrals of motion). Additional physical requirements can be used 
to fix solutions. For example, in the case of a perfect fluid the flux A is fixed at the critical 
point [21[3]. 

Integrating Eqs. (fT4"j) and (fT5|) we find, 

rr 

M(V, r) = C + AV + 4tt / T °(r)r 2 dr, (17) 

Jr 

where Co is a constants of integration and r$ is the horizon locus. Similarly, integrating 
Eq. (fTTI) . in the limit A € 1, one obtain 

/•T 

A = -47T / T^rdr + siV), (18) 

J r 
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s(V) is an arbitrary function. Note that integration of Eq. (|12p gives the same result for 
A, Eq. (|18p . In Eq. (|18p the function s(V) can be set to zero by redefinition of the time 
coordinate, e s ^ v 'dV — > dV. The constant Co, in turn, can be identified with Mq by an 
appropriate shift V — > y+const. Therefore Eqs. ()1T|) and f)18|) can be rewritten as follows, 

rr 

M(V,r) = M + AV + 4tt T °(r)r 2 dr, (19) 

J ro 

rr 

A(r) = -4vr / T x °rdr. (20) 



Jr 

The above equations are the central result of this work. 

If the components of the energy-momentum tensor are slowly varying functions of the 
radial coordinate (which is a reasonable assumption for non-pathological matters), then 
from ()19p and (|20p one can find the expressions for the corrections of the metric in the 
vicinity of the black hole horizon, 

M(V,r) = M + AV + 4irr$ (r-r )T ° \ r=ro , (21) 
A(r) = -4irr (r - r ) \ r=ro . (22) 



According to our perturbative scheme, the above equations (|19p and (|20p are only valid 
if the corrections are small. Thus we require, 



\AV\ < M , 



,r 

4tt / T °(r)r 2 dr 

J ro 



< M , 



I- / T^rdr 



< 1. (23) 



'r 

The first of the inequalities in (I23p is simply a condition that the total flux of the infalling 
matter is smaller than the bare black hole mass. Similarly, the second condition says that 
the total mass of the matter inside a radius r is smaller than the bare mass of the black 
hole. The last condition in (|23p has no obvious interpretation, it reflect the fact that the 
correction to the non-diagonal part of the metric due to the presence of infalling matter 
must not be too big. 

2.3 Reissner-Nordstrom black hole 

We can generalize our results to the case of the Reissner-Nordstrom black hole. Accord- 
ingly, instead of ([8]) we now write the metric coefficient v in the following form, 



where Q is a charge of a black hole. We assume that the accreting fluid carries no electric 
charge, otherwise the total charge of the black hole needs to be modified as well. Sub- 
stituting (j24"|) into (f9|)- (fT2|) and taking into account a contribution of the electromagnetic 
field in the total energy- momentum tensor, 

rp . rp , Q 2 rp 1 rp 1 . Q 2 



8vrr 4 ' 1 1 87rr 4 ' 

V^T--^, r ,._ >T ,"-^ l , (25) 

after lengthy but straightforward calculations one can check that the expressions for the 
metric coefficients M(V, r) and X(V, r) are the same, as in the case of the Schwarzschild 
black hole, and given by Eqs. (|19p and (|20p . where ro has now the meaning of the event 
horizon of a bare charged black hole. 
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2.4 Accretion of a perfect fluid 

In this section we consider accretion of a perfect fluid onto a black hole. The stress tensor 
of a perfect fluid is given by, 

Tftu = (p + p)u^Uu ~ P9^u, (26) 

where p is the energy density and p is the pressure of the fluid in the comoving coordinates, 
and is a 4-velocity of a fluid. In the Eddington-Finkelstein coordinates the four-fluid 
velocity is given by 

u" = ( _L , -u, 0, ] , (27) 

where u = \dr/ds\ > is the absolute value of the radial component of 4-velocity in the 
static coordinates and /o = 1 — 2Mo/r. Note that the value of u is finite at r = tq. It is 
also easy to check that all the components of vP" and are non-divergent at the horizon. 
The relevant components of the energy-momentum tensor are 



PV fo + u 2 ~pu „ o P + P 



r u = Ti u = (28) 



In the vicinity of the horizon, /q — > 0, we find, 

T °^1(P-P), T,^- P -±f. (29) 
Thus, from (|21|) and (|22p we can write for the corrections, 

Af(V,r) w M + AV + 2vrrg (p - p) (r - r ) , (30) 
A(r) « ^ ro ^( r _ ro ). (31) 

The above expressions are valid for any perfect fluid in the vicinity of the horizon of a 
black hole. Note that the energy flux onto the black hole is given by (|16p . with Tq 1 = 
— (p + p)u\f /o + u 2 for a perfect fluid, see [3]. Therefore, as it can be seen from the 
(f30|) . the accretion of phantom energy, p + p < 0, decreases the black hole mass. Thus 
we confirm that taking into account backreaction does not affect the results, which were 
reported previously in [3], where the zeroth-order approximation has been used. 

Note that the expression (130p implies, in particular, that light observers, having dif- 
ferent V = const, see the black hole with the different masses M (at different stage of 
accretion). 

2.5 Shift of the apparent horizon 

The position of the apparent horizon depends on a choice of a coordinate system. For the 
metric ([6]) it can be shown that the location of the apparent horizon, r^, can be found 
from the following equation [25], 

e viy ' r) = 0. (32) 

The above expression can be obtained by requiring that for a radial outgoing photon 
dr/dV = 0. Indeed, from ds 2 = one can find the two radial light geodesies, 

dr 1 

dv = < W e '"" = r < 33 > 
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The apparent horizon satisfies the condition that the photons do not cross the r = const 
surfaces in the growing r-direction. It gives dr/dV = 0, and from ()33[) one obtains e^ +A = 
and hence e v ^' r ^ = 0. The later is due to A regularity. Now, for a Schwarzschild black 
hole from (f32|) and (JT9J) we find, 

M + AV + A-K [ k T °(r)r 2 dr = ^, (34) 
Jro ^ 

which is an implicit equation for r^. It is not difficult to check that for small shifts of the 
apparent horizon the last term in the l.h.s. of (|34p can be neglected, since it is of the next 
order in the expansion, so that we find, 

r h « 2M + 2AV. (35) 

Thus the leading term in the apparent horizon shift only depends on the total flux A and 
is independent on other components of the energy-momentum tensor. 

Similarly, for the shift of the apparent horizon in the case of a charged black hole, one 
has, instead of (fM)) . 

r h ru o 2 

M +AV + Att T °{r)r 2 dr = ^- + ^-. (36) 

Again, neglecting the last term in the l.h.s. of the above relation, we obtain, for the case 
of the Reissner-Nordstrom black hole, 



r h » M + AV + yjM 2 -Q 2 + 2M AV. (37) 

Note that for the extreme black hole, Mo = Q, the existence of the apparent horizon 
depends on the sign of the flux, A. For phantom accretion, the apparent horizon ceases 
to exist for positive V, which means that our perturbation scheme breaks down for this 
particular case. The reason for the breaking down of our scheme is that accretion of 
phantom decreases the black hole mass, therefore any amount of phantom matter converts 
the extreme black hole to a naked singularity, thus going beyond our quasi steady-state 
approximation. On the other hand, if normal (non-phantom) matter accretes, then (|37p 
can be perfectly applied. It is interesting to mention, that in contrast to the Eddington- 
Finkelstein coordinates, which we use here, in the static Schwarzschild coordinates the 
test fluid approximation is violated for accretion of any type (phantom and non-phantom) 
matter onto the extreme black hole 121 1. 



3 Relation to other solutions 

Having found the general form of the corrections to the metric due to backreaction, we 
can now compare our results to other solutions. In particular, it is interesting to see how 
our result is related to known full analytic solutions. 

Let us first consider the case when matter is the cosmological constant term. In this 
case it is easy to verify that T ° = px, 3\ = and A = 0, where pa is the energy density 
of the vacuum term. The general expression (|19p and (|20p give then, 

M(V,r) = M + ^-p A {r 3 -r 3 ), A = 0, (38) 
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which, after the redefinition Mq — > Mo+47r/?A^o/3, can be identified with the Schwarzschild- 
de Sitter solution. 

As a second example, we consider the case of radially infalling photons, so that the 
energy-momentum tensor is of the form, 

V = -1^2^- ( 39 ) 

Then the only non-zero component of the energy-momentum tensor which contributes to 
the corrections of the metric (|19p and (|20p is Tj corresponding to a non-zero total flux, 
A^O, while T ° = T 1 ° = 0. Thus, 

M(V, r) = M + AV, A = 0, (40) 



which corresponds to the Vaidya solution [16] for the constant flux, A =const. 

It is also interesting to observe that in the vicinity of the horizon of a Schwarzschild 
black hole, independently of the form of the energy-momentum tensor of the accreting 
fluid, the corrections to the metric are of the Vaidya form. Indeed, for r = 2Mq, from 
CEHD and flU, we infer that M(V,r = r ) = AV, A = 0, which coincides with (00]). This 
means, in particular, that the first correction to the motion of free particles in the vicinity 
of the horizon does not depend on the energy-momentum tensor of the accreting fluid: it 
depends on the total flux A only. One should note, however, that although the metric 
itself is asymptotically Vaidya near the horizon, the derivatives of the metric are not the 
same (not even asymptotically) as in the Vaidya solution. Said differently, the Riemann 
tensor is different — even near the black hole horizon — for different accreting matter. 

4 Examples 

4.1 Stiff perfect fluid 

In the zeroth-order approximation steady-state accretion in the Schwarzschild coordinates 
corresponds to a steady-state solution in the Eddington-Finkelstein coordinates. Therefore 
we obtain the same integrals of motion (as a functions of r only), the same expression for 
the critical points, and finally the same profile of density and velocity as in both coordinate 
systems. Thus we can take results of [3], derived in the static coordinate system and apply 
them here. In particular, for the linear equation of state p = p — Po one has 

,=£ + ( fa «u 1+ ?w 1+ 3y (4i) 



4 V 4 / V v J V 7*^ 

where poo is the energy density at the infinity. The corresponding profile for the velocity 
u = u(r) can be calculated from the *** energy flux *** conservation equation. In 
the case of a stiff perfect fluid (po = 0), we find, see [3], pu 2 = roPoo/V 4 . The relevant 
components of the energy-momentum tensor in this case are given by, 

0_„ (, , r 0\ A r o\ r r0 _ ,„ ^ 



T u = poo(l + ^(l-^j, ^ = -2^(1+'^). (42) 

and the flux, 



A = 8n Poo rl (43) 
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Therefore the corrections in the vicinity of the horizon take the form 



M = Mq + %TT Poo r$V + 47rp 00 rg I— + — - x 2 + - - log xj , (44) 
A = iirpoorQ (x 2 + Ax - 5 + 2 log x) , (45) 
where we have introduced a short-hand notation x = r/rp. 

4.2 Scalar field 

Here we consider accretion of a scalar field and find corrections to the metric. The bound- 
ary conditions at the infinity are posed by the cosmological evolution, 

<A|oo = 4>c- (46) 

Let us first consider the standard kinetic term (which in the hydrodynamical description 
correspond to a stiff perfect fluid), so that the action is given by 

Scan = / d*Xy/=jj (\d^(ft) (47) 



2 

Solving the Klein-Gordon equation, □</> = 0, in the metic of the Schwarzschild black hole 
one can deduce that the stationary non-singular solution for <ft is given by, 



r-r-rolog(^) ). (18) 



Taking into account that the energy-momentum tensor for a canonical scalar field is 

T il y = d il <ftd v <ft- 1 -{d<ftf 9tiv , (49) 

one can easily calculate the relevant components of the energy momentum tensor on the 
solution (gS]), 

V = + V = -^(l + 2!) 2 . (50) 



2\ r I \ r J \ r 



A = 4rr$rg. (51) 



and 



Substituting the above expressions to (119]) and ([20]) one can see that the correction to the 
metric are given by (|44p and (|4"5]) after identification (ft 2 . — > 2p OQ ^ as one could expect. 

It is worth to mention phantom accretion in the context of backreation effects. One 
of examples for phantom, which was considered in |3J, consisted of a scalar field with the 
"wrong" sign in the action. In terms of a perfect fluid this simply means that the density 
and the pressure of the effective perfect fluid becomes negative. This, in particular, implies, 
from (|44p that the black hole mass decreases during accretion of a phantom field. This 
result takes into account backreaction of the accreting field and it confirms the conclusion 
of |3|, where the backreaction was neglected. 

Now let us consider accretion of a galileon scalar field. The general galileon model 
consists of five non-trivial terms: the first one is a potential of the form const Xcft, the 
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second one is the canonical scalar field (which we considered already above) and three 
other non-linear terms. Here we consider accretion of the galileon with the action, 



al 



The stress-tensor for this action is, 

iff = -24>,^ V U4> + 2<t>^V u) {d<pf - g^ a V a (d^y 



(52) 



(53) 



The solution for steady-state spherically symmetric accretion onto a black hole was found 
in 13, 

' r 



(f) = 4> c \ V - r + - 2^/r^log 

Substituting (f54"|) into (f53"j) and ([TUP one finds for the galileon, 
; o 3 „n : o 3 



+ 1 



: x(l + ^)' 1 
Thus the correction to the metric are 

M = Mo + Utt^V 

+ Airfcrl ( 6 (v^ - 1) - 3 (x - 1) + 2 



(54) 



(55) 



.3/2 



1) -61og(v/x + l 



A 



127T^r 2 



1 + 



+ 2^-41og(l + v / x)-l + 41og2) . 



(56) 
(57) 



5 Conclusion 

In the studies of accretion of matter onto black holes, there are basically two analytic 
approaches. One is to neglect backreaction of an accreting matter onto the metric of the 
black hole — which is reasonable when the matter is light — and to think of the matter 
as a test fluid on top of the background metric. Another way is to try to find solutions for 
the full system of the Einstein and matter equations. There are only a few such solutions 
known, the reason for this is that it is in principle very hard to find solutions to the full 
non- linear system of equations. 

In this work we studied the backreaction of an accreting fluid on the black hole met- 
ric in a perturbative way. Accretion is spherically symmetric and quasi steady-state. In 
this approximation we assumed that the accretion flow is small and the accreting fluid is 
light, so we could use the perturbation scheme. At zeroth-order approximation the accret- 
ing matter does not backreact, so this step is equivalent to the test-fluid approximation, 
Eqs. ([3]) and (j3J). Then the first-order corrections to the metric are found from the per- 
turbed Einstein equations, in which the source term is given by the accreting matter in the 
zeroth-order approximation, Eq. ([5|). In particular, for the case of spherically symmetric 
accretion onto a Schwarzschild black hole, the resulting equations for the correction are 
written down in ([9l)- (ll3p . The solution for the corrections to the metric coefficients are 
found in general form, Eq. (|19p and (|20p . this is the central result of our paper. In a similar 
manner, the general perturbation scheme is applied to the Reissner-Nordstrom black hole, 
and we have found the same result for the correction of the metric, Eq. (|19p and (|20p . 
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We compared our perturbative results with known analytic solutions. In particular, we 
verified that the our scheme gives correct results in the case of the cosmological vacuum 
term (up to a redefinition of the central mass), as well as for the radially infalling photons, 
in this case we reproduced the Vaidya solution. 

An interesting consequence of the general result is universality of the corrections to the 
metric in the vicinity of the horizon of a black hole. We found that, if the accretion rate 
is nonzero, then, independently on the energy-momentum tensor of an accreting fluid, the 
leading corrections to the metric is of the Vaidya type, Eq. (|40p . 

We applied the obtained results for calculation of correction to the metric for a generic 
perfect fluids and also for particular examples: accretion of a stiff perfect fluid Eqs. (|44p. 
(I45p as well as accretion of a canonical scalar field and the galileon scalar field, Eqs (|56p . 

dSZD. 
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